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Abstract: We analyze the axially-symmetric scalar perturbations of 6D chiral gauged super- 
gravity compactified on the general warped geometries in the presence of two source branes. 
We find all of the conical geometries are marginally stable for normalizable perturbations (in 
disagreement with some recent calculations) and the nonconical for regular perturbations, 
even though none of them are super symmetric (apart from the trivial Salam-Sezgin solution, 
for which there are no source branes). The marginal direction is the one whose presence is 
required by the classical scaling property of the field equations, and all other modes have 
positive squared mass. In the special case of the conical solutions, including (but not re- 
stricted to) the unwarped 'rugby-ball' solutions, we find closed-form expressions for the mode 
functions in terms of Legendre and Hyper geometric functions. In so doing we show how to 
match the asymptotic near-brane form for the solution to the physics of the source branes, 
and thereby how to physically interpret perturbations which can be singular at the brane 
positions. 
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1. Introduction 



Six dimensional supergravity continues to provide an instructive laboratory for exploring 
different aspects of the cosmological constant problem 0, |2|. It does so by building on the 
earlier observation ||| |], [5|, ||] that higher-dimensional brane constructions can break the 
implication that says that a nonzero 4D energy distribution must necessarily generate a large 
4D curvature. Once recast into a higher-dimensional context the usual cosmological constant 
problems j/J become questions as to what are the choices required for obtaining acceptably 
flat 4D space, and identifying how natural these choices may be. Part of the appeal of the 6D 
proposal is that it is very predictive, having observational implications both for accelerator 
physics and for tests of gravity in the lab and in astrophysics Q , in addition to its potential 
implications for late-time cosmology Q. 

In extra-dimensional models there are two separate kinds of naturalness issues which 
must be addressed. The first type of naturalness demands stability under renormalization 
for any choices which are required to ensure the flatness of 4D space. That is, any choice 
of brane positions or couplings which are required for 4D flatness must be stable against 
integrating out the ultraviolet part of the theory. Although extra-dimensional models appear 
to make progress in understanding this kind of stability |lCj ] , it remains open whether they 
can completely resolve this part of the problem. 

The second type of naturalness issue arises in extra-dimensional models because these 
models generically relate the Dark Energy density to the dynamics of the various moduli 
fields which describe the geometry of the internal dimensions. Although it is possible to 
obtain phenomenologically successful cosmology from this kind of dynamics ||, the second 
naturalness issue asks whether this requires unusually delicate choices for the brane and bulk 
initial conditions. In particular one worries that generic initial conditions — or small per- 
turbations to specially-chosen ones — might give rise to catastrophically fast time evolution, 
leading to late-time cosmology which is unlike what we see around us. 

The present paper addresses the second of these naturalness issues. Any understanding 
of this issue requires knowledge of the time-dependent solutions of the system. And since the 
time scales of interest are not slow compared to those associated with the Kaluza Klein states, 
this time evolution is intrinsically higher-dimensional and cannot be accurately described 
within an effective 4D description Instead, what is required is an intrinsically higher- 

dimensional description of the system's dynamics and its response to various probes. 

In a companion paper |^] we present a broad class of time-dependent scaling solutions 
for 4D compactifications of 6D supergravity. These solutions aid in understanding the evo- 
lution produced by generic brane properties by showing what the late-time evolution of the 
system is likely to be for a wide class of initial conditions. The present paper is aimed at 
a complementary part of the problem of time-dependence. Instead of seeking the late-time 
behaviour towards which the system evolves, we here study how small perturbations can 
initiate the beginnings of time-dependent behaviour. (See ref. for a description of some 
of the intermediate and transient time-dependent phenomena which lies between these two 
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cases.) We do so by investigating the linear perturbations of a broad class of static solutions 
to 6D supergravity, in order to see how they respond to arbitrary small perturbations. In this 
analysis, we use a Kaluza Klein decomposition which includes both the modes oscillating in 
time as well as the zero-mass non-oscillating ones. 

Our main result is that these systems are marginally stable for normalizable perturbations 
around conical geometries, and for regular perturbations around nonconical ones. (We also 
examine a class of physically motivated non-regular perturbations, for which our stability 
analysis is inconclusive.) The stability is only marginal because of the generic scaling property 
which the 6D supergravity solutions have. (It is the absence of an energy barrier in this 
marginal direction which allows the development of the scaling solutions studied in ref. fl2|| .) 
We identify the set of equations which govern the time evolution of the system after all of 
the constraints have been used to eliminate some of the field fluctuations, and reproduce the 



form found by earlier workers [14]. But we differ in the conclusions which we draw regarding 
stability from these equations, for reasons we explain in detail in the Appendix. 

We begin our discussion in §|| with a reminder of the field equations and some of the 
known static background configurations. This is followed in §|| by the derivation of the 
equations governing linearized perturbations in comoving gauge, while the same derivation in 
longitudinal gauge is performed in appendix [A| for comparison with previous results. then 
solves analytically the linearized equations in the conical case and identifies the asymptotic 
behaviour in the more general non-conical case. We then give in §[| several general arguments 
in favour of the stability of the modes examined and relate them with the boundary conditions 
in §||. After a brief summary of conclusions - §^ - four appendices contain: details of the 
analysis in longitudinal gauge, with the full stability argument; a detailed comparison of why 
our stability conclusions differ from some previous results; the relation between some of the 
results in various gauges; and a summary of some properties of the special functions used in 
the explicit solutions. 

2. Field Equations and Background Solutions 

We begin in this section by summarizing the relevant field equations, as well as the static 
background solutions about which we shall perturb. 

2.1 Field equations 

The bosonic part of the Lagrangian density for 6D chiral gauged supergravity is given (for 



the case of vanishing hyperscalars — < 5 a = 0) by [15 1 1 

2g 2 



C 1 

_ g MN 



-g 2k 2 



Rmn + d M cpd N (p 



K 4 



\ e-* F MN F MN - -L e-^ G MNP G MNP , (2.1) 



lr rhe curvature conventions used here are those of Weinberg's book |16[ , for which all curvature tensors 
differ by an overall sign relative to those of MTW ITS . 
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where F MN = d M A N — d N A M and G MNP = d M B NP + d N B PM + d P B MN + (A P F MN terms). The 
coupling constants g and k respectively have dimension (mass) -1 and (mass) -2 . 
The field equations obtained from this action are: 



b 



2g 2 



D M [e- 2 * G MNP 
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+ d M <pd N <p + ^ e" 2 ^ G„poG w pq + K 2 e~* F ArP F, p 
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+ 2 (° ¥>) ^ 







0. 



(dilaton) 
(2-Form) 
(Maxwell) 

(Einstein) 
(2.2) 



An important feature of these equations is their classical scaling property. This property 
states that given any solution to these equations another can be obtained by making the 
replacement 

e v^\e>f and g MN ^ X'^mn , (2.3) 
with all other fields held fixed. This property follows from the fact that the action, eq. (|2.1|), 



scales under this transformation according to S — ► \~ 2 S [18]. 

The remainder of the paper deals with solutions to these equations for which two di- 
mensions are compact and four are not. To this end denote the six coordinates x M , M = 
0,1,2,3,4,5, as x M = {t,x,y,z,i],6}. We denote 4D coordinates by x M , \x = 0,1,2,3, 
where x^ = {t,x,y, z}, and 2D coordinates of the extra dimensions by x m , m = 4,5, where 
x m = {r/,0}. When required, the three spatial coordinates of the observable four dimensions 
are denoted x l , i = 1, 2, 3 and so x l = {x, y, z}. 

2.2 Ansatze 

In later sections our interest is in metrics of the form 

d s 2 = e 2a rj^ dx»dx v + e 2v dr, 2 + e 2b d6 2 (2.4) 

and consider these component functions, as well as (p and Ag, to depend only on the coordi- 
nates rj and t. (In what follows we sometimes generalize this assumption to allow dependence 
on rj and x^.) Denoting differentiation with respect to rj and t by primes and dots, the field 
equations for these functions then reduce to the following set of coupled partial differential 
equations. 

The Maxwell equation is: 



„2(v-a) 



A e + (2a + v-b- y)A e + A' e ' + (4a' - v' - b' - <p')A' e = . 



(2.5) 



The dilaton equation is: 



a 2(v-a) 



cp + (2d + v + b)ip + if" + (4a' - v' + b')tp' 



1 

+ 2 6 



e -2a+2v-2b-<p /j^}2 + ± p ~2b- v( A l^2 _ p 2v+ V 



(2.6) 
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The (trj) Einstein equation is 

3d' + 1/ - 3va' + bbf - bd - vb' + ifxp' + e -26- * 'A e A' e = 0. 
The (tt) Einstein equation is: 



(2.7) 



D 2(v-a) 



3 a + b + v + {vf + (6) 2 + (ip) 2 - a(v + b) + a" + a' (4 a' - v' + b') 



3k 2 



-2a+2«-26-^ )2 _ «_ ^26-^)2 + 5 e 2,+^ _ Q (2 _ 



The (r/r/) Einstein equation is: 



D 2{v-a) 



i\2 



v + v{2a + ij + b) + 4 a" + 6" + 4 (a') 2 + (&') - t/(4o' + 6') + (^) 



-2a+2v-2fc-<^ 



+ 



— e - 2b -^(A' e ) 2 + ^e 2v+ v = 0. 



(2.9) 



The (00) Einstein equation is: 



O 2{v-o) 



b + b{2a + v + b) + b" + b'{4a' -v' + b') 



3k 2 



-2a+2v-2b-f i 



(i,)2 + e -™-v(A' e ) 2 + 4 e 2v+ v = . (2.10) 



The (ij) Einstein equation is: 



2(v-a) 



a + a(2a + v + b) + a" + a (4a - v + b') 



+ ^ e -2a+2v-2b- v ^ g ) 2 - !j e- 2b ^(A' e ) 2 + ^ e 2v+ ^ = . (2.11) 
2.3 Maximally symmetric 4D compactifications 

Many explicit compactifications of the above field equations to four dimensions have been 
constructed over the years, starting 20 years ago with the Salam-Sezgin spherical solution [15]. 
These now include compactifications to flat 4D space on unwarped, rugby-ball solutions [Q], 
as well as warped axially-symmetric internal dimensions having conical Q, and more general 
[ |Te| , 20, 21] singularities at the positions of two source branes. More recent generalizations 
have found warped de Sitter 4D geometries [22], as well as configurations for which the 
hyperscalars and 3-form fluxes are nontrivial p3[ |. 

Since we wish to perturb about the 4D flat solutions, we briefly summarize their properties 
here. 2 These solutions may be written 



AW 4 



and 



e b = A, 



2 The conventions of ref. may be obtained from ours by taking 7?a 
= 1/2, while those of §) differ from those here only by the choice k =1. 



-Rmn, f 



(2.12) 

-y/2 and 
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and so the bulk fields become 



ds 2 = W 2 (ri) ri^v dx^dx" + A 2 (t]) >V 8 (r?) dry 2 + d# 2 
Frfi = e~ x ^ and e~f = W 2 e x ^ , 



(2.13) 



where 



d 



k 2 <?A 2 \ cosh[Ai(r/ - r/i)] 



2gAi / cosh [A 2 (77 - m)] 

• 2„3 



,4 " 4 = (l^f ) e " 2A3 ' ?cosh3 [ A i( r ?- ? ?i)]cosh[A 2 (7 ? - ??2 )], (2.14) 

where the field equations imply the constraint A 2 = A 2 + A 2 . 

There are two particularly interesting special cases of these solutions. The first is obtained 
by taking A3 = 0, in which case the resulting geometries have purely conical singularities [2C]. 
The second comes by taking both A3 = and 7/1 =7/2, in which case the conical geometry 
is also unwarped, since <p = <po and W = Wo are constants. In this latter case changing 
variables to proper distance, dp = AW 4 dr], the Maxwell field becomes 

F pd = d p A e = ±%e^B(p), (2.15) 

where B(p) = e b<yp ^ satisfies 

^ = 6" + (6') 2 = -^e^, (2.16) 

with solution B(p) = Bosm(p/po) with po = ^fte~ v °/ 2 / g. These are the unwarped, rugby- 
ball generalizations Q of the original spherical Salam-Sezgin solution jll| to include the 
back-reaction of two branes, located at the sphere's north and south poles. 

In what follows we also use the extremely useful change of variables used in ref. [22], 



<p = l(X-y-2Z), logW = j{y-X) and logA = ^(3X + y + 2Z), (2.17) 



and so 



e -x = e-ip/2^-1 = cosh [ Al ^ _ m y (2.I8) 
Ai 

25 



e~y = e-^A^W" 4 = -f- cosh [A 2 (r? - m )} (2.19) 

kA 2 



e 



- 2 = e v W 2 = e~ x ^ . (2.20) 



It is straightforward to check that these are related by the Hamiltonian constraint (for evo- 
lution in the 77 direction): 

e ^y + K y e 2* + x 12 - y 12 + Z' 2 = 0. (2.21) 

K 
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2.4 Asymptotic forms 

The solutions to these equations describe geometries which typically become singular within 
the bulk. Because of the assumed axial symmetry this can occur at two (or fewer) points, 
which we can choose to place at r\ = ±00. These are interpreted as being the positions of 
the 3-branes which source the corresponding field configuration. We require an explicit form 
for the asymptotic behaviour near these singularities in order to determine the boundary 
conditions in our later linearized perturbation analysis. We therefore pause here to outline 
what this asymptotic behaviour is. 

If we use proper distance as the radial coordinate near the branes we have ds 2 = 
g a b dx a dx b + dp 2 , and we imagine the brane position to be given by p = 0. For maximally- 
symmetric 4D geometries the asymptotic form of the bulk fields in the p — > limit is then 



generically given by a power law form [22] 



ds 2 ~ [c a {Hp) a f q^ dx^dx" + dp 2 + [ce(Hpf) 2 d9 2 

~ c^Hpf and FP e ~c f (H P y, (2.22) 

where a, (3, p, 7, c a , c$, and c/ are constants, H is an arbitrary dimensionful parameter, 
and q^y denotes the metric of 4D flat, de Sitter or anti-de Sitter space. For time-independent 
fields the field equations impose the following two Kasner-like constraints amongst the powers 
a, j3, 7 and p: 

4a 2 + p 2 +p 2 = Aa + (3= 1 and 1 = p -l. (2.23) 

Imposing these 3 conditions on the 4 powers leaves 1 undetermined, which we can choose to 
be p. Together with the prefactor, eg, this power can be related to the two physical choices 



which are available for the source branes: their tension and dilaton coupling [T^, 20, 24, 25]. 
For example, explicit calculation using the solutions given above gives 

A2-A1 A 2 + 3A lT 2A 3 , 2(A 2 -A!T2A 3 ) 

a± = — : : — , P± = — : — and p± = — , (2.24) 

5A 2 - Ai t 2A 3 5A 2 - Ai t 2A 3 5A 2 - Ai T 2A 3 V ; 

for the asymptotic form as r] — > ±00. As is easily verified, these satisfy eqs. ( |2.23| ) given 
above. 

For time-dependent geometries it is useful to generalize the above to allow different 
asymptotics for the spatial and temporal components of the metric, replacing eq. ( 2.22j ) with: 



ds 2 ~ -[c w {H P T] 2 dt 2 + [c a (Hp) a ] 2 g^dM + dp 2 + [c e (Hpf] 2 d6 2 , (2.25) 

with the powers now related by 

u? + 3 a 2 + p 2 +p 2 = uj + 3a + p = 1 and j = p - 1 . (2.26) 

The additional parameter appearing in this asymptotic form is related to the freedom to 
separately specify the energy density and pressure on the source branes. 
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3. Linearization 



In this section we set up the equations which govern axially-symmetric perturbations about 
the solutions described above which transform as scalars in 4D. The restriction to axially 
symmetric perturbations does not limit the ensuing stability analysis because the most un- 
stable modes are generally the most symmetric, since any angular dependence contributes 
positively to the corresponding Kaluza-Klein mass. 

3.1 Symmetries and gauge choices 

The most general 4D scalar perturbations have the form 

ds 2 = e 2a [r,^ + M, dx^dx" + N m ^ dx»dx m + g mn dx m dx n 
B = B^ dx^ A dx v + B m>li dx^ A dx m + B mn dx m A dx n 
and A = Q >fi dx M + A m dx m , (3.1) 

where in four dimensions B^ v dualizes to a 4D scalar £ through a relation of the form H^ u \ oc 
e^\ p d p (. We are free to use gauge symmetries to locally set N m = B m = £2 = 0. 

Discrete symmetry and mode mixing 

As has been pointed out by earlier analyses ]14|, p6fl , for the purposes of a stability analysis 
it is not necessary to keep all of the remaining perturbations. To see why, notice that since 
our interest is in fluctuations depending only on t, x l &nd r/, the absence of a dependence on 
the coordinate 9 implies also a symmetry under the reflection, 9 — ► —9. If Ag is nonzero 
in the background configuration (as is the case when a flux F v g is turned on in the extra 
dimensions), this is only a symmetry if we also independently require A M — > —A M . 

This symmetry ensures that the linearized fluctuations can be divided into two classes, 
according to whether or not they are even or odd under these reflections. In particular we 
have 

(even) 



{5a, 5ggg,5g m ,5<p, 5A e } 

{Sgrfi^SB^SAJ (odd). (3.2) 



Since the symmetry guarantees that it is consistent with all of the equations of motion to 
set the odd fluctuations to vanish, these two categories of fluctuations cannot mix with one 
another at the linearized level. We take advantage of this fact to choose bg^Q = b(, = SB^g = 
SAr) = 0. Again, the omission of these modes does not restrict the stability analysis which 
follows. 

We are led in this way to the following ansatz for the metric and other bulk fields 

ds 2 = e 2a (r,^ + M M „) dx»dx u + e 2v di] 2 + e 2b d9 2 
A = A e d9 and £ = . (3.3) 
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To linearize the field equations about a static background we write 

a(rj,x) = ao(rj) + A(r),x) , v(rj,x) = v (rj) + V(rj,x) , b(rj,x) = b (rj) + B(rj,x) , 

(f(r},x) = (p (ri) + $(r/,x) and A e (r), x) = a e (rj) + A e {i], x) , (3.4) 

where we allow the fluctuations to depend on all 4 coordinates, x^, and the background field 
configurations are given as above: 

e ao = W(rj) , e v ° = A(v)^(v) , e b ° = A{rj) , 

e -w = W 2 (??) e x sn and fl / = e -*v . (3.5) 

For some of the discussion to follow it is useful to choose proper distance computed using the 
background metric as a coordinate within the extra dimensions, as we did for the rugby ball 
solution above. We reserve the variable p for proper distance, and it is related explicitly to rj by 
dp = e v °^dri = AW^dn. Notice that this choice of coordinates can be made independently 
of the gauge choice we make on the linearized fluctuations, which we now describe. 

Three useful gauge choices 

Finally, it is also not necessary to keep all of A, M, V and B independent, because one 
combination of these can be set to zero using an appropriate coordinate choice. There are 
three gauge choices which we use in what follows. 

• Comoving (c) gauge: defined by the condition A^f 1 = 0, leading to 

ds 2 = e 2ao e 2A( C ) + e ^o e 2V^ ^2 + ^2 + ^ _ ( g g) 

together with A = agd6 and cp = ipo + <&( c \ Notice that this gauge does not fall 
completely into the ansatz of eq. ( |2.4| ), because the coordinate transformation required 
to get into this form does not preserve the condition 5Ag = 0. This is the origin of the 
new metric variable, fow, in the above expression. 

• Gaussian Normal (GN) gauge, defined by V^ GN ^ = and so: 

ds 2 = e 2ao (e 2A(G "V + Mf v N) ) dx^dx" + e 2b0 e 2BiGN) de 2 + dp 2 , (3.7) 

together with ip = ip + $( GAr ) and Aq = a$ + A g GN \ It is obviously convenient also to 
use background-metric proper distance, p defined via dp = e v °dr]. 

• Longitudinal (I) gauge: defined by M = and so 

ds 2 = e^e^n^dxW + e 2b ° e 2B<l> d9 2 + e 2y( °dp 2 , (3.8) 

together with (p = <pq + $W and Aq = oq + Ag . Notice that we again choose to use 
background-metric proper distance, p, as the radial coordinate. 
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Which of these gauges is most convenient depends on the goal of the calculation. Co- 
moving gauge is physically intuitive (since one follows perturbations along hypersurfaces of 
constant Ag), and has the enormous benefit that the field equations completely decouple in 
this gauge, making the linearization analysis much easier. Longitudinal gauge is convenient 
for making contact with earlier calculations in the literature, since this is the gauge which 
has often been used. Finally, GN gauge is useful for analyzing the boundary conditions which 
fluctuations must satisfy near the source branes. 

It is clearly useful to be able to transform between expressions obtained in these three 
gauges, so we next briefly summarize the main formulae which are required. To change from 
longitudinal to comoving gauge requires the transformation = xfy + £ M , with £ M = e5™ 

so that pa\ = pr c \ + e, with e = A® /d p ag required to ensure A^ = 0. This implies that the 
various fields in longitudinal gauge are related to the ones in comoving gauge by: 

A( l ) = + e d p a , = + e d p b , F« = + d p e , (3.9) 

AK C ) = _ e 4ao+6o £ ; $(0 = $(c) + edpipo 

Similarly, to go from comoving to GN gauge, we perform the transformations p( c ) = 
P(GN) + e w flh dp £ = chosen to ensure V^ GN ^ = 0. The further transformation, = 
x^Q^+d^e, with d p e = e~( 6a ° +b °^ — e~ 2a °e is also required in order to ensure the vanishing 
of N( GN \ Quantities in comoving gauge are then related to those in GN gauge by 

A (c) = A (GN) + ^ ) M (GN) = _ 2£) B (c) = B (GN) + ^ 

$ (c) = QiGN) + £dp ^ A f N) = -ed p a e . (3.10) 
3.2 Linearized equations in comoving gauge 

We now work out the perturbed equations in comoving gauge, since in this gauge we are able 
to decouple the fluctuations even when perturbing about non-conical background solutions. 
We use as the coordinate r] in this section, with the derivative d/dry denoted by primes. 

We adopt the following notation for the perturbations A^ = — \l//2 and = £/2, and 

so 

ds 2 = eP-^e-^d^d^+e^^ ( 3. n) 

while e v = e (x-y-2Z)/2 e ® an d p^g = a' e = qe~ X3T >e 2X+z . In terms of these variables, and 
Fourier transforming with respect to time and space 3 so □ = Tf v d p d u = uj 2 , the Maxwell 
equation is 

2 uj 2 e {x - y ^ 2 N + 2B' + 2$' + £' + 4*' = ; (3.12) 



3 We emphasize that ui is the mass of the corresponding Kaluza-Klein mode and should be understood as 
the eigenvalue of the operator □. In this analysis, we do not make the assumption that the modes have a 
specific time-oscillatory behavior of the form e lfe °*. Thus this analysis includes both time-oscillating modes, 
and non-oscillating ones which correspond to zero-mass Kaluza Klein modes. 
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the dilaton equation is 

<j>" + ^^+2$ - I L 2 q 2 e 2X + ^e 2y ^ $ = <f/ + ^' ~ & + wV*-^/ 2 ^ 

+ ^ e ^ + K 2^ e 2AT s (313) 

the (nrj) Einstein equation is 

2{2X' + Z')B - {2y' + Z')i + + AB' - 6*' = , (3.14) 

while the other off-diagonal term in the Einstein equation (the combination (ii) — (jj), with 
i / j) is 

e f(#-5y-2Z) ^iv' + (Af' - y)JV] - 25 - £ + 2^ = 0. (3.15) 

A combination of the (7/77) Einstein equation, the dilaton equation of motion and the trace of 
the Einstein equation which involves only first derivatives with respect to 7/ is 

^ e (*-y)/2 (2 y + 2 > )N _ V e 2y^ + 3^ e 2y+zy + (2 y + ^ 
2 k 2 

The (00) Einstein equation is similarly 

(25" + $") + uj 2 e 2y+z {2B + *) - 2n 2 q 2 e 2X {2B + $) + 2A"#' 

= A" + £' + 2 a; 2 e ( * _y)/2 iv) , (3.17) 

and finally, the (/J/x) Einstein equation is 

($" + ^ e 2y+z <$>) - (V + u 2 e 2y+z * - 2Z'¥) = Z' (b' - ±? - M*- y V 2 N) .(3.18) 

We have here a set of seven equations for five variables (£, N, B and <£). Two of these 
equations appear as constraints and can be used to fix the lapse and shift functions: using the 
(/xry) equation ( 3.14] ) and the (7777) equation ( |3.16| ), we can express both £ and u 2 N in terms 
of the three remaining variables ^f,B and <E>, leading to five remaining equations for three 
variables. However these remaining equations are not independent since the Bianchi identities 
ensure that two of them are redundant. One can be expressed as a linear combination of the 
others, whilst the derivative of the Maxwell equation ( |3.12 ) is simply the derivative of ( |3.14 ) 
and a combination of the other equations. 

We arrive in this way with three equations for three unknowns, which can be written as 

1 r,^ , , ^, , 16 g 2 e 2y 



B " + u 2 e 2y+z B + \Ux' + 2y> + 4Z> + B> 
4 \ k 2y' + Z' J 



+\ [~l2. 2 q 2 e 2X + ^g^) B + \ (3*' + y + 2Z>) # (3.19) 
1 ( 10 „2J^ , !63 2 o y 2X'-3Z'\^ 6g 2 e 2y 



._ -\2K 2 q 2 e 2X + e 2y — — $ - ^ - ^' = 

16 I q k 2 2y + Z> 2Y + Z' ' 
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* + „V>«* + ^ p^*' + \tf - Y) (*< + 2BQ - ^y^B' (3.20) 

K 2 2 o 2 p 2y 

e 2X (* + 25) - 1JL [(2*' - 3Z')* + 2(2A" + = , 

and 

4 n 2 p 2 y 

$" + w 2 e 2 ^d> + y/ (*' + 25') + I [W - 2B') (3.21) 

K 2 n 2 2 n 2 p 2 y 

e 2X (* + 25) - f*. [(2Af - 320* + 2(2Af + = . 

These equations dramatically simplify once the change of variables (<£, 5, — > (f, x, $ 
is performed, with 

Kq 

V 1 2U 2 + X'Z' x ~ , 

B = -— x --* + e ~ x f 3.22 

2 2nqU 

T 2y + 7 2^ 2 - 3yz' _ - A"y ~ 

4 V 4V3WV 2 Kg 6/ 

where use of the background field equations shows that the nominally field-dependent quan- 
tities IA and V are really both constants: 

U 2 = X' 2 + K 2 q 2 e 2X = X 2 

V 2 = X' 2 + K 2 q 2 e 2X + \z' 2 = X 2 + |a§ . (3.23) 

Notice that for conical backgrounds IA = V = Ai. In terms of these new variables the three 
field equations — eqs. (|3.19|) , ( |3.20|) and ( |3.21|) — take the remarkably simple, decoupled 
form 

x" + » 2 e 2 y+ z x = (3.24) 
/" + uj 2 e 2y+z f - U 2 f = (3.25) 
i2o2y+z ~ _ 3V 

(2y + z'y 



+ w 2 e 2 *+^ - °-±JL- -JLl _^ = o . (3.26) 



Because these equations decouple from one another each can be solved separately, and 
presents no problem for numerical analysis. Exact solutions are also available in the case 
of backgrounds having conical singularities (including but not restricted to the rugby ball 
geometries), as shall be described in detail in later sections. Given the existence of variables for 
which the equations decouple in comoving gauge, one might wonder what the corresponding 
variables are in other gauges. These variables can be obtained by explicitly performing 
the required change of gauge, but because the results are not very simple we relegate their 
description to Appendix |C[ 
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It is noteworthy for later purposes that the decoupled equations ( p. 24 - 3.26 ) are inde- 
pendent of the background variable X (and hence also of the parameter rji). Consequently 
the linearized equations of motion for the new variables x, f and tp are identical for a general 
conical background geometry (for which Z' = 0) and a rugby ball geometry (for which Z' = 
and X' = y'). However once the solutions are inserted into the expressions for the metric 
fluctuations using relations Q3.22 ), the results do depend on X and so the full fluctuations 
'know' about the complete background geometry. 



4. Properties of Solutions 

Up to this point there has been no restriction on the form of the background geometry. In this 
section we specialize to special conical backgrounds for which it is possible to solve the above 
linearized equations in closed form. We then return to the more general non-conical geometry 
in order to discuss the asymptotic behaviour of the solutions in the near-brane limit. The 
section closes with an outline of how this asymptotic form is related to the physical properties 
of the branes which source the bulk geometries. 



4.1 Analytic Solutions for General conical backgrounds 

We now concentrate on the behaviour of the perturbations to solutions having only conical 
singularities, which are defined by the condition A3 = (and so Z' = 0). Performing the 
change of variable z = tanh [A2 (ij — ^2)], the perturbation equations become 



fx 

dz 2 
d 2 j> 
dz 2 

dz 2 



2z 



2z 



dx 
dz 

d$ 
d7 



+ 



+ 



/i- 



2z 



1 



d/ 

Tz + 





-z 2 


(r 


1 




-z 2 


(r 






-z 2 



i 1 



4i'' 



1 



f = o, 



(4.1) 
(4.2) 
(4.3) 



with /j, 2 = uj 2 k 2 /Ag 2 , and to is the 4D mode mass. Notice that these equations are symmetric 
under z — * —z, which corresponds to reflections of 77 around 7/ = 7/2- The reason for this 
symmetry is clear in the special case of the rugby-ball solutions (for which rji = 772), since 
it then corresponds to reflections of the spherical geometry about its equator. In the more 
general warped conical geometries this symmetry instead follows from the above-mentioned 
circumstance that the equations are completely independent of r/i , and so take the same form 
for general conical geometries as they do for the rugby ball. 

-1 + \J\ + 4/i 2 — or equivalently writing /x 2 = z/(i/ + l) — the solutions 



Defining v 



to these equations become 
X = C 1 P u (z) + C 2 ReQ u (z) 



= C X F 



■ z 



„ ^I> + 1) 1 



(4.4) 

1 / s v 3 9 
-(l + u),l + -;u + -;z- 2 
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/ 1 



1/ X V 



1 



vT 



CrF 



2' 2' 
1 



+ C 4 z 2 F 



2'"2 (1 + ^ ); 2 ;Z 



+ C 6 z F 



v 1 
2' 2 
1 



v 3 
2' 2 : 



(4.5) 
(4.6) 



where the Cj, i = 1, ..,6 are integration constants and F(a,b; c; z) is the hypergeometric 
function. Neglecting the overall normalization of all three modes, we see that the above 
expressions provide a 9-parameter family of solutions. Notice also that these solutions can be 
decomposed into symmetric and anti-symmetric combinations under the symmetry z — > —z, 
and this is already explicit for the functions ip and /. 
We see that the mode energies satisfy 



(4.7) 



with 1/ro = 2<7/k, and so the spectrum depends on the values which are allowed for the 
parameter v. Notice that for real v we have u> 2 > except for the interval — 1 < v < 0. As 
usual, the values which are allowed for v are related to the behaviour which is demanded of 
the solutions at the boundaries, which in the present instance corresponds to the near-brane 
limits z — ► ±1. For this reason we next explore the asymptotic form taken by the solutions 
in the near-brane limit. 



4.2 Asymptotic Forms 

We next examine three aspects of the near-brane limit. In this section we first explore the 
asymptotic forms taken by the solutions just constructed for fluctuations about backgrounds 
having only conical singularities. We then see how these asymptotic forms generalize to 
the broader situation where the background has more generic singularities. A preliminary 
connection is drawn in the next section between this asymptotic behaviour and the properties 
of the source branes. 



Conical Backgrounds 

The asymptotic form of the solutions found above near the brane positions at z = rbl may 
be found by using the properties of the Hypergeometric functions given in appendix [D|. This 
shows that for generic z/, x an d ip diverge logarithmically as one approaches the branes, 
whereas / generically diverges as {l-z 2 )- 1 / 2 ^ cosh[A 2 (7? -772)] 

As usual, less singular asymptotic behaviour can be possible for those v for which the 
hypergeometric series terminates, which for F(a, b; c; z) occurs when either a or b vanishes. 
Inspection of the explicit solutions then shows that termination happens when v is quantized 
to be a non-negative integer: v = £ = 0, 1,2, . . .. Notice that the resulting quantization of 
KK mass, from eq. (4/7), implies in this case u> 2 > 0. 

For v = I one - but not both - of the two Hypergeometric functions appearing for each 
of x, ip and / terminates, allowing a less singular near-brane limit to be obtained through an 
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appropriate choice for the integration constants Cj. For instance, if we require — as we shall 
argue on the grounds of normalizability below — that all fluctuations must be less singular 
than (1 — z 2 )" 1 , then this is possible for / only if v = I and one of C5 or Cq vanishes (which 
one depends on whether £ is even or odd). 4 Once this has been done the functions x an d 
ip can still diverge logarithmically as z — > ±1, leaving a 4-parameter family of potentially 
logarithmically singular solutions (up to overall normalization). 

Completely regular solutions are also possible when v = £, provided that one of C\ or C2 
and one of C3 or C4 also vanishes. For instance, x ls regular everywhere \iv = £ and C2 = 
since in this case the solution Pi degenerates to the Legendre polynomials. This leaves a 
2-parameter family of nonsingular solutions (up to overall normalization). 

Non-conical backgrounds 

Although it is difficult to find explicit solutions in the generic case where the background 
has non-conical singularities. It is possible to use the linearized equations to estimate the 
asymptotic behaviour of the mode functions as ij -> ±00 (or z — > ±1). Repeating the steps 
which led to the equations of the previous section for this more general case, we instead find 



fx 

dz 2 

d 2 j> 
'dz 2 

fl 

dz 2 



2z 



l-z u 
2z 



1 



dz (1 



dip 



n 2 



z 2 ) \1 



l + z 



6\ 



dz 



+ 



2z 



df 

Tz + 



A 2 



x = o 

'l + z 



(\-z 2 
,2 



(1 



1 - 2 
l + z 



1 



dA 



6\ 



2 (1 - SX 2 
(z - 5\) 2 

1 - 45 A 2 



/ = 



(4.8) 
$ = (4.9) 

(4.10) 



with fl 2 = /i 2 e x ^ 2 and SX = A 3 /2A 2 . 

Since \5X\ < 1/2 the term in eq. ( [4.8D which is proportional to x remains subdominant 
compared to the terms in (l — z 2 ) 1 dx/dz and d 2 x/dz 2 , such that the asymptotic behaviour 
of x still remains logarithmic near the brane. Similar arguments show that ip and / also keep 
the same asymptotic behaviour as in the conical case. 

Normalizability 

A restriction on the asymptotic forms follows from the condition that the mode functions be 
normalizable, as we now explore. As we show in detail in the next section, the equations of 
motion for each of the decoupled modes can be derived from a reduced action, eq. ( |5.11|) . In 
this action all modes, u, have the same kinetic term, e 2 ^ +z uDu, where e 2 ^ +z = .A 2 W 6 , and 
so the definition of the inner product can be taken to be the same for all three of Xi V' an d /, 
having the form: 

(u,v) = - [ d 3 xdedrjA 2 W 6 {ud t v* -u*d t v), 



4 We give explicit expressions for some of the lowest of these nonsingular modes for ip and / in Appendix A 
for longitudinal gauge. 
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for any two modes u and v. Here S denotes a surface of constant t. The norm of a mode is 
therefore defined by 



(u,u)=uj f d 3 xd6dr/A 2 W 6 uu* = u [ d 3 x d6 dpAW 2 uu*, (4.11) 

where we use dp/drj = AW 4 and dtu = —iu>u. 

Our interest is in the conditions placed on the asymptotic behaviour of linearized modes 
by the convergence of these integrals. Given the asymptotic form, eq. ( p. 22 ), of the back- 
ground metric, with powers ao, Po, Po etc., we have AW 2 oc pP° +2a< -i ^ pi-2a j n ^he ne ar- 
brane limit (p — > 0). In order to be normalizable, modes should therefore vary as -u oc p v with 
power t; > — 1 + ao as P ~~ * 0. In particular, for a conical background with Qo = and Po = 1, 
modes must diverge less rapidly than 1/p near both branes. 

4.3 Boundary conditions and brane properties 

The above considerations show that the stability of the linearized fluctuations is related to 
the behaviour of the modes as they approach the branes. In this section we briefly explore 
to what extent the properties of the linearized solutions capture what we know about the 
asymptotic properties of the exact solutions. Exploring this connection also allows us to 
make some contact between the conditions for stability and the physical properties of the 
branes which source the bulk field configurations, through considerations along the lines of 



those given in refs. [12, 2C, 24] 



Asymptotics in GN gauge 

Recall for these purposes that on general grounds the various bulk fields are known (in GN 
gauge) to approach the branes through a power-law form given by eqs. ( |2.22| ). In principle 
there are two cases to consider, depending on whether or not the background and perturbed 
fields share the same powers in these asymptotic expressions. That is, if we imagine that the 
background and perturbed metric components, e a ° and e a = e a °e A , satisfy 

e a ° - c a0 (Hp) a ° and e a ^ c a (Hp) a , (4.12) 

near p = 0, then we see that the fluctuation field, A, must satisfy 

A -»• (a - a ) ln(Hp) + In f — ) + • • • , (4.13) 

in the same limit. Similar expressions also hold for the fields B, V, and 

Since the Maxwell field strength varies as F p q oc pT+ 2 ^ oc pP +2 ^~ 1 the same argument has 
slightly different implications for Aq, since clq oc 

pPo+2{3 and A e = a e +A d (x pP +213 . (Notice 
that this gives the expected results F p $ oc p and Aq oc p 2 in the conical case, for which (3 = 1 
and p = 0.) In this case we have Ag oc pP +2 P if p + 2/3 < po + 2(3o, because the small-p 
behaviour of Aq dominates that of ag. Things are different if p + 2(3 > po + 2/3q, however, 
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since in this case ag dominates Ag , and in order to achieve this Ag must contain a piece which 
varies in the same way and cancels the contribution of ag as p — > 0. As a consequence in this 
case Ag oc pP°+ 2 ^ for small p. 

How singular this is depends on the allowed range of p + 2(3. Notice that in general the 
constraints 4a + (3 = 4a 2 + (3 2 + p 2 = 1 require \(3 2 — \(3 + p 2 = \- Solving this equation for 
p, the sum p + 2/3 is thus minimal (' — ' sign) and maximal ('+' sign) for f3 = g =F ' an( ^ 
p = =F-^, so that the sum is in the range -1.433 ~ §(1-V2l) <P+2(3< §(l+>/2l) ~ 2.233. 
We see that perturbations to .Ag can become quite singular. In the particular case of conical 
singularities p = and (3 = 1 and so p + 2/3 = 2, leading to smooth behaviour near the branes. 

We see from the above that the perturbations need not be smooth at the brane positions, 
but that in GN gauge the metric perturbations can at worst diverge logarithmically in the 
near-brane limit. (Perturbations to the Maxwell field can be more singular than logarithmic, 



but only if p and (3 are sufficiently negative.) Since the arguments of refs. (T^, |2C], 24] relate 
powers like a to physical properties on the branes, these logarithmically singular perturbations 
only arise if some of the brane properties are themselves perturbed. It is only in the limit 
that the asymptotic near-brane behaviour is the same before and after perturbation that the 
perturbed solutions are nonsingular at the brane positions. This is the case in particular if it 
is assumed that the geometry has purely conical singularities before and after perturbation, 
or if the sources are represented as delta functions with the implicit associated assumption 
that the bulk fields are well-defined when evaluated at the brane positions - as is common in 



higher-codimension stability analyses, such as that of ref. [14|. 



Asymptotic forms in comoving gauge 

Although we have seen that within GN gauge fluctuations diverge at most logarithmically 
near the source branes, we next determine what this implies for the near-brane behaviour in 
the (comoving) gauge we use in our analysis. In this section we examine in particular the 
near-brane asymptotics in comoving gauge. 

Choosing background proper distance, p, (rather than rj) as the radial coordinate we have 
seen that the perturbations in comoving and GN gauge are related by eqs. ( |3.10| ), which we 
reproduce here for ease of reference: 

A (c) = A (GN) + £dpaQ j B (c) = B (GN) + £dpbQ 

V {c) = d p e and $ (c) = $ (G7V) + ed p ip , (4.14) 

with e = —Ag GN ^/(d p ag). The condition = is ensured through a change of the 

coordinates with parameter e = —M^ GN ^ /2. 

Since we know how the fluctuations vary with p as p — > within GN gauge, the above 
equations allow this information to be carried over to comoving gauge. In particular, since 
ag + Ag oc pP +2 P and ag oc fP 0+2 "° , we saw Ag oc where Q = min(p + 2(3, po + 2/?o). We see 
that e oc p 1+A and that there are two cases: (i) A = (p — po) + 2{(3 — (3q) if p + 2(3 < po + 2(3$, 
or (ii) A = if p + 2(3 > po + 2(3o. Since ao, bo, (po, vary at most 
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logarithmically for small p, we see from the above that all of the fluctuations in comoving 
gauge are also at worst logarithmic provided only that A = (jp + 2(3) — (p + 2/3 ) > 0. 
(In particular, for conical backgrounds the perturbations are at most logarithmic provided 
A = (;p + 2/?-2)>0.) 

These arguments show that the fluctuations <£, B and \& diverge at most logarithmically 
within comoving gauge for a broad range of physical situations. Inspection of the definitions, 
eqs. (3.24)-(3.26), shows that this implies the same conclusion for e~ x /, x and 'tp. 



5. General Stability Analysis 

We now give several general arguments in favour of stability for a broad class of boundary 
conditions. We provide two types of arguments, which complement one another. The first of 
these works directly with the linearized equations derived above, and can be used directly with 
either the longitudinal-gauge (provided in Appendix (El)) - or comoving-gauge formulations. 
The second argument is instead cast in terms of the action, and closes a loophole left by 
the previous equation-of-motion analysis. We make this second argument only for comoving 
gauge (and not for longitudinal gauge, say) due to complications which arise in constructing 
the relevant action in other gauges. 



5.1 Equation of motion and tachyons 

We start by arguing for stability directly with the equations of motion. The goal of this 
argument is to relate the sign of the energy eigenvalue, u 2 , to the boundary conditions which 
the fluctuations satisfy near the positions of the source branes. 

The most direct way to do so is to multiply eq. ( |3.24| ) by x*; sum the result with its 
complex conjugate; and integrate the answer over the extra dimensions, to get: 



to 



/oo 1 r i oo r+oo 

^e 2y+Z \x\ 2 = -\\{\x\ 2 )'} +/ drylxf. (5.1 
-oo Z L J -oo J_ oc 





oo f 


\\x\ 2 r 


+ / 




-oo J_ 



This shows how the sign of lu 2 is related to the behaviour of \ near the two brane positions. 
Notice that the combination e 2y+z appearing on the left-hand side of ( |5.1[ ) is related to A 
and W by e 2y+z = A 2 W 6 . 

Although we do not have explicit solutions to the linearized equations for general non- 
conical backgrounds, the asymptotic behaviour of these solutions was argued above to be the 
same as in the conical case. In particular, both fluctuations x an d ip are typically proportional 
to logp ~ T) as rj — > ±oo, whereas / usually diverges as 1/p ~ in this limit. Thus the 
general asymptotic behaviour as rj — * ±oo is 

X~ (C±r] + D ± ) , ^~ (E±r] + F±) , 

and / ~ (G±e X2H + H+e'^A , (5.2) 
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for constants C± , D± , E± , G± and H± . Using this in eq. (|5.1|) , and cutting off the integrations 
at r\ = ±A, gives 



j +K drie 2y+z \x\ 2 = -{\C + \ 2 + \C-\ 2 )A-Re(D* + C + -D*_C^)+ f + d V \x'\ 2 ■ (5.3) 
J-A J-A 



We see from this that if the fluctuation is required to remain finite on both branes — i. e. C+ = 
C- =0 — the boundary term vanishes and the squared energy of any mode is necessarily 
positive: to 2 > 0. Furthermore, lo only vanishes if x is a constant throughout the entire bulk. 
Applying the same argument to the -0 equation gives: 



co 2 I A d7 ? e^M 2 = -- (|,f)' +/ d// Mr? ' ^ \r\ 
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32 g 2 V 2 e 2y 



k 2 (2y + z 



(\E + \ 2 + \E_\ 2 ) A-Re(F|£+ - F1EJ) (5.4) 

[ +K A A 7/|2 32 g 2 V 2 e 2y . 
+ / dr? |V / | 2 + J , 



A 



(2y + Z') 2 



Here again, u 2 > if £/+ = E— = is imposed so that the fluctuation remains finite on both 
branes, and can only vanish if ip vanishes identically. 

Finally, the same argument applied to the mode / gives 

- 2 £ A dr/e^l/l 2 = -\[(\f\ 2 )'} + _ A A + jT A dr? (|/'| 2 + U 2 |/| 2 

= -A 2 (|G+| 2 + |G_| 2 ) e 2A2A + A 2 (|# + | 2 + \H^\ 2 ) e~ 2A2A (5.5) 

so that finiteness of the boundary terms requires G± = if A 2 > or H± = if A 2 < 0. In 
either case lo 2 > and can only vanish if / = throughout the bulk. 

In this way we see that perturbations about any background (let it be conical or not) 
are marginally stable, provided that the fluctuations are restricted to be nonsingular at the 
brane positions. 

5.2 Action analysis and ghosts 

The previous argument shows that the eigenmode frequency lo must be real for a broad choice 
of boundary conditions at the brane positions. However one might worry that the Lagrangian 
density for the linearized fluctuation might have the form C = $*/((/?)[□ — g(op)]&, where <p 
(resp. $) denotes a generic background (resp. fluctuation) field, and f(ip) and g((p) are 
background- field dependent quantities. Notice that the equation of motion for $ implies 
[— □ + g((p)]$> = 0, and so implies lo 2 > if g((f) > 0. But this might still imply a negative 
contribution to the fluctuation energy if the prefactor f{ip) should happen to be negative for 



some configurations ip. (Indeed, precisely this form of instability was argued in ref. [14] to 



occur in 6D chiral supergravity linearized about rugby-ball configurations due to kinetic-term 
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mixing amongst various KK modes.) We next present a second argument for stability, which 
closes this particular loophole. 

The starting point in this approach is to expand the action, eq. (|2,1|), to second order in 
the perturbations. Working in comoving gauge and using the coordinate rj we get in this way 
the following quadratic action 
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•^kin 



J d 6 x (£ kin + £ dyn ) 
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4k 2 



r e2 
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where we allow the fluctuations to depend on all four noncompact coordinates, x^, □ = 
Vj^d^dy denotes the flat d'Alembertian, and the split into 'dyn' and 'kin' distinguishes terms 
which involve □ from those which do not. We obtain the above expression by freely integrat- 
ing by parts while ignoring the surface terms which this procedure introduces at the brane 
positions. This neglect of surface terms is justified inasmuch as our goal is to exclude the 
possibility of the negative (ghost-like) kinetic terms, as described above. 

Differentiating this action with respect to the shift function TV" leads to the constraint 

2{2X' + Z')B + 4<p' $ + 2 (2B' - 3*') 



^ 2y + Z' 

while differentiating with respect to the lapse function £ gives the constraint for iV: 
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These two constraints are consistent with what was obtained in eqs. (glj) and ( gig) . As 
before, we may use both these constraints to obtain the equations of motion for the three 
remaining degrees of freedom, and Vl/. Substituting eqs. ( |5.9D and ( 5.10| ) into the action, 
together with the change of variable ( |3.22 ) , the resulting action takes the remarkably simple 
form 
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(5.11) 
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Varying this action gives the same equations as the ones obtained in ( 3.24 , 3.25 , 3,26 ), confirm- 
ing that the substitution of ansatz ( |3.6| ) into the action consistently reproduces the equations 
of motion of the full theory. 

Since all of the kinetic terms in the action are positive, it follows immediately that the 
theory has no ghost modes (for which f(<p) < 0). Because this disagrees with the conclusions 
drawn in [14], we also reproduce their longitudinal-gauge analysis in Appendix [B] and identify 
a sign error which we believe to be the source of the discrepancy. 



6. Conclusions 

In this paper we have studied the linearized evolution of perturbations to gauged chiral 6D 
supergravity compactified to 4D on a broad class of static and axially-symmetric vacuum 
solutions, sourced by two space- filling 3-branes. We follow previous authors in focussing on 
scalar modes which share the axial symmetry of the background, and which are even under 
a convenient parity transformation. These restrictions are consistent with the equations 
of motion, since they are enforced by symmetries. They do not compromise the stability 
conclusions because the excluded modes necessarily have higher squared-frequencies than do 
the lowest of the modes kept (which we find are bounded below by uj 2 = 0). 

Although none of the backgrounds we perturb are supersymmetric, we find they are 
all marginally stable to perturbations which are well-behaved at the brane positions. The 
marginal direction is unique and is the one required on general grounds by a general scaling 
property of the 6D supergravity equations. Besides providing a general stability argument 
which works for a broad class of static vacua, we also provide analytic solutions to the fluc- 
tuation equation for the special case of conical backgrounds, including (but not restricted 
to) the rugby ball, (in comoving gauge) and the rugby ball solutions (in longitudinal gauge). 
These allow us to explicitly verify that the near-brane behaviour of the solutions has the 
properties required by general asymptotic properties of the field equations. By constructing 
the truncated action for the lowest-lying KK modes, we are able to show that no instabilities 
arise from mode-mixing amongst the lowest-lying levels, contrary to recent claims. 
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A. Analysis in longitudinal gauge 



In this section, we linearize the field equations in longitudinal gauge, specializing to those 
background configurations for which A3 = and so which only have conical singularities. 
Working in this gauge allows us to make contact with earlier calculations, and so for ease of 
comparison in this section we adopt the notation of these earlier workers. Following refs. 26, 
|L"f | we therefore write 



A (i) 



"2 ' 



2 ' 



$(0 



— and Jip 



(A.l) 



We also follow refs. |^6|, Q and adopt as radial coordinate, p, as proper distance in the extra 
dimensions (measured using the background metric), leading to the full ansatz 



ds 2 = W 2 (p) e"* v„u dx»dx v + e« dp 2 + A 2 e 2 *^ d9 2 



(A.2) 



together with e~ 2lfi = W 4 e* and A' e = (qA/W 6 ) + A e . Throughout this appendix we use 
primes to denote d/dp rather than d/drj. 

Fourier transforming in time allows us to write —d 2 = to 2 , and so we find that the (tt) 
and (ij) Einstein equations degenerate into the same equation at the linearized level, to give 
the following {pis) Einstein equation: 
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Similarly, the (tp) equation becomes: 
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The (pp) Einstein equation becomes 

AW' 2 A' 



W A\ , 



2W 



f-2K 2 



2W 2 2 ^ 

3k 2 



,2 



w 



10 



A 

2* - 



& + - 
2 



1 fAW 3A 



W 



+ 



A 



2 2 



AW 4 



A, 



AW 4 



A 6 



2W 



9 



W 

2 



2 W 2 



f 



0. 

(A.4) 

(A.5) 
= 0. 



The (99) Einstein equation becomes 
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(A.6) 
= 0. 
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The dilaton equation becomes 

u 2 f .„ (AW\ ( Tl j\ /4W A! 



and the Maxwell equation becomes 



w 2 y V w .4/ v w6 7 V 2 

Nominally we have here 6 equations — dilaton, Maxwell, and the (tt), (tp), (pp) and 
Einstein equations — for the 4 unknown functions, ^, £, / and ^4#- However, as was 
the case in comoving gauge, the Bianchi identities ensure that two of these equations are 
not independent. This is because two combinations of these equations can be chosen to only 
depend on single derivatives with respect to time, and the Bianchi identities ensure that the 
solutions to these lower-order constraint equations is consistent with the time evolution. This 
allows us to use the constraint equations to solve for two of the fields, and then to drop two 
of the other field equations as redundant. 

Following this logic we determine Aq by using the (tp) Einstein equation, which integrates 

to 



V-4W 4 



\a 1 fr> ^ ( W A '\ t ( W A '\, ( W \ r 



where Cl(p) is an arbitrary function of integration. This unknown function can be determined 
by using the equation obtained by adding the (tt) + (pp) + (86) Einstein equations (which is 
a constraint inasmuch as it does not contain any time derivatives): 

Substituting eq. ( A. 9 ) into eq. ( A.lOj ) and using the background equations of motion leads to 
the following condition for 0: 

(AW A n)' = 0, (A.ll) 



and using this in the derivative of eq. (|A.9|) gives 

-2„ \ i /qvw' a'\ /3W' 3.4'V, (W 
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We may now use this last expression to eliminate A'g from the remaining 3 field equations, 
which we take to be the dilaton and (tt) and (pp) Einstein equations, leading to: 
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and 



W 2 ' s \w / \W 

/6W' „4'\ 8 5 2 / A 

in agreement with eqs. (37)— (39) of ref. [jl4|]. These equations easily lend themselves to 
numerical integration because they are independent, since all of the constraints have been 
explicitly solved. 

The above equations take a somewhat simpler form if rewritten in terms of new variables: 
x = 2f + / = -4A - 2$ = 2(B + V - $) and T = V + f = 2(V - A) = B + 3V: 

+ x + hn- + T x=°. ( A - 16 ) 



W 2 V w 



and 
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w 2 V w 4 7 w 2 V ^ 



The logic to solving these is to integrate eq. (A.16|) for \i an d then to use the result as a 



source in eqs. ( A. 17] ) and ( |A.18 ) to integrate these. These equations may be solved in closed 
form for the case of the rugby-ball solutions, as we show in what follows. 

Although we arrive at equations which are in perfect agreement with ref. p4|], we shall 
differ in the conclusion we draw from them inasmuch as we find the system to be marginally 
stable (modulo an issue concerning boundary conditions which is irrelevant for the comparison 
with pjl). For this reason we also review the stability analysis of ref. jl4j in appendix ^ and 



explain the origin of the discrepancy. 

A.l Analytic Solutions for the Rugby Ball 

The rugby ball is the special conical geometry for which e w ° = e a ° = W = e^ = 1. In this 
case explicit forms for the axially-symmetric scalar mode functions can be found analytically 
in longitudinal gauge in terms of well-known special functions. This has the advantage of 
providing a concrete verification of the previous stability arguments, as well as providing a 
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framework within which to see how the boundary conditions give rise to the quantization of 
Kaluza-Klein frequencies. 

Recalling that in this section, we use the proper distance, p, as radial coordinate the 
rugby-ball geometry is given by 

e bo = A = A r sin ( — ) where r = — . (A.19) 
V r o/ 2 9 

These choices imply A = 1 — 6 encodes the defect angle, 2ir5, at the branes situated at p = 
and p = tttq. 

Using the rugby-ball conditions W = and A' /A = Vq 1 cot(p/Vo) allows the fluctuation 
equations in longitudinal gauge, eqs. ( |A.16| )-( |A.18| ), to be simplified to 

-W' + ^^x' = o 



w , 3cot(p/rn) 2T 
u 2 T + T" + ^— ^ r' g- 
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(A.20) 



2 f , f // cot(p/r ) _ 2 2cot(p/r ) , 
uj J -\- J j — — w l x 



To solve these equations we write uj 2 = p? ' jr\ and change variables to z = cos(p/ro), leading 
to a relatively familiar set of equations 

d 2 X ( 2z \ d X ( p? \ 
dz 2 \l-z 2 ) dz^~ \l-z 2 J X 

^ - ^ + r - -JL_ ( A 2i) 

dz^ I 1 - Z 2 J dz + { ! _ Z 2 ) 1 - ! _ Z 2 



dz 2 \1 — z 2 J \1 — z 2 J \1 — z 2 J dz 

In terms of the variable z = cos(/o/ro) the brane singularities correspond to z = dbl, at 
which points we have 

z = 1 — — -) , near p = 

z = -1+ fc~^ r °) 2 + ... , near p = vrr , (A.22) 

and so p — p± = =F\/2(1 ± 2:) 1//2 + • • • where p+ = 7rro and /9_ = 0. The asymptotic near-brane 
boundary conditions are therefore related to the behaviour of the solutions as z — * ±1. The 
asymptotic form which is appropriate to longitudinal gauge may be found by performing the 
transformation from GN gauge, along the lines as was done in the main text for comoving 
gauge. In this section we concentrate on finding those solutions for which all perturbations 
are less singular than (1 — z 2 )^ 1 oc p~ 2 , including in particular those for which the divergence 
at z — ► ±1 is at most logarithmic. 
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Homogeneous solutions 

We start by providing explicit integrals of the homogeneous parts of eqs. (A.21), for which 
the right-hand sides vanish. The first of these is the Legendre equation and therefore has as 
solutions 

X (z) =C l P u (z) + C 2 ReQ u (z), (A.23) 

where v{y + 1) = fi 2 , or v = v± = |[— 1 ± (1 + 4^/ 2 ) 1//2 ]. We take the real part of Q u here 
because this function is normally complex for \z\ < 1. (See Appendix ^ for our detailed 
conventions.) Since V- = — 1 — v+ and the differential equation is invariant under v — > — 1 — v 
it suffices in what follows to use only the positive root, v = ^[-^/l + Ap? — 1]. 

For generic v both P v and Q u have logarithmic singularities as z —* ±1. However, we 
shall find that for generic u, T diverges like (1 ± z) . Restricting to at most logarithmic 
solutions therefore requires us to take v = I = 0, 1, 2, . . ., and so we only consider these values 
from this point on. (However, we emphasize that since logarithmic singularities are physi- 
cally acceptable for the present problem, this quantization is not required by the boundary 
conditions for X-) The mass formula 
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(A.24) 



then shows that stability directly follows from this quantization, in agreement with the general 
arguments of earlier sections. 

Given v = I the hyper geometric series for Pi{z) terminates and degenerates to the Leg- 
endre polynomials, which are regular at both z = ±1. Nevertheless, since Q v {z) has the 
asymptotic limit 
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(A.25) 



it clearly retains its logarithmic behaviour at z — ► ±1 even when v = t. Combining the 
above expressions allows the near-brane singularities to be written in terms of the remaining 
integration constants, C\ and C 2 , as follows: 



X{z) 



C2 
2 

C 2 



ln(l -z) +0(1) 



as z 



1 



2(-l) 



In (1 + z) + O(l) asz^-1. 



(A.26) 



The homogeneous part of the equation for V is also of Hypergeometric form, leading to 
the following general homogeneous solutions: 
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Notice that the function multiplied by C3 (or by C4) is even (or odd) under reflections about 
the 'equator' at p = vrrg/2, since these correspond to z — > —z. In the special case v = the 
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Hypergeometric series can be summed explicitly to give 



F[l,b;b;z 2 ]=J2' 
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(A.28) 



k=0 



showing that in this case the two homogeneous solutions are the elementary functions (1 
z 2 )- 1 and z^l-z 2 )- 1 . 

The asymptotic behaviour of the Hypergeometric functions as z — > ±1 turn out to be 
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which shows in particular that diverges like (1 



(A.29) 

2 ) _1 for generic v. Since we seek 
solutions which only diverge logarithmically, we must choose either C3 or C4 to vanish, and 
then choose v so that the Hypergeometric series for the other solution must terminate. That 
is, either C4 = and \ — ^ = — n (and so v = 1 + 2n) or C3 = and 1 — § = — n (and so 
1/ = 2(n + 1)), where n = 0, 1, 2, ... is a non-negative integer. These are automatically satisfied 
for non-negative integers, v = £ 7^ 0, with even (odd) £ corresponding to C3 = (C4 = 0). 
For v = we have seen that the homogeneous solution is (C3 + C4 z)/(l — z 2 ), which is only 
nonsingular at both z = 1 and z = — 1 if C3 = C4 = 0. We see in this way why it is that v 
must be quantized to be a non-negative integer. 

The nonsingular homogeneous solutions for for the first few choices for £ are: 
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£ = 1 implies = C3F 
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and so on. 

The homogeneous part of the equation for / is also Hypergeometric, and is given in this 
case by 



f h (z) = C 5 (l-z 2 )F 



+ 2'2 2'2' 
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Notice that these are the same Hypergeometric functions as appear in T, and so their singular 
properties can be read off from the expressions given above for T^. In particular, because of 
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the additional factor of (1 — z 2 ) the asymptotic forms do not diverge as z 
Hypergeometric series does not terminate, leading to the limits 
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For the first few choices for I this leads to the following homogeneous solutions 



(A.31) 
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and so on. We see that because of the quantization, v = £, the only logarithmic singularities 
in the solutions are those appearing through Qe(z) in \- 

Perturbing tensions 

Dividing the perturbations into those which are even or odd under the reflection z — * —z allows 
an even more precise interpretation. Finite perturbations which are even under this reflection 
correspond to perturbations which change the tension of both boundary branes in the same 
way, and so remain within the class of rugby-ball solutions (for which both brane tensions 
must be equal). Conversely, regular perturbations which are odd under z — > —z correspond 
to perturbations which change the two brane tensions oppositely, and so describe excursions 
into the more general class of warped, conical solutions [l9|, Logarithmically diverging 
perturbations describe the transition from the conical to nonconical class of solutions. 



Particular integrals 

We now return to the problem of solving the differential equations for T and /, including 
now the right-hand sides of these equations. Since the general solution is found by adding 
the general solution to the homogeneous equation to any particular solution, it suffices here 
to identify particular integrals to these differential equations, which we denote by T p i(z) and 
fpi(z). The general solutions are then simply T(z) = T^z) + T p i{z) and f{z) = fh{z) + fpi(z). 
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For instance, if we specialize to nonsingular solutions at both z = ±1 we may choose 
C2 = and so take x to be given by the lowest few Legendre polynomials. This leads to the 
following expressions for T(z) = Th(z) + T p i(z): 
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(A.32) 
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and so on. 

With these solutions in hand we next find the relevant particular integrals for the differ- 
ential equation for /. Again restricting to nonsingular solutions, specializing to the first few 
Legendre polynomials, and summing f = fh + fpi gives: 
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C 4 z; 



and so on for as many modes as are desired. 

A similar process can be undergone in the more general case where C2 7^ 0, and so 
X = C\ Pi(z)+C2 Re Qi(z). In this case x diverges logarithmically in the near-brane limit, and 
the particular integrals are no longer obtainable in such a simple closed form. Nonetheless they 
may be obtained numerically in terms of integrals over the right-hand-side of the equations 
weighted by appropriate combinations of Hyper geometric functions. 

A. 2 Stability Analysis 

Using the explicit previous solutions and particularly their asymptotic behaviour, we give 
here a general argument in favour of stability for a broad class of boundary conditions. 



As performed in section 5.1, we argue for stability directly with the equations of mo- 
tion and relate the sign of the energy eigenvalue, u> 2 , to the boundary conditions which the 
fluctuations satisfy near the positions of the source branes. 

Proceeding as before, we multiply eq. ([A. 16 ) by AW 4 x*'i sum the result with its complex 
conjugate; and integrate the answer over the extra dimensions, to get: 
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We denote here the brane positions by p = and p = p\. Since \f—g = AW 4 > we see 
that lu 2 > follows provided we can show that x satisfies boundary conditions for which the 
combination AW 4 d(|x| 2 )/dp vanishes at the brane positions. 

In general the appropriate boundary conditions in longitudinal gauge for x have x ~^ 
(Hp) Al in the near-brane limit, p — > 0, with A; being a power which is determined in terms 
of the asymptotics of the background configuration, and A; = could correspond to pertur- 
bations which remain finite at the branes as well as those which diverge only logarithmically 
at the brane positions. Using this asymptotic limit we have 
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(A.35) 



dp dp 

as p —* 0. This shows that the integral on the left-hand side generically diverges due to its 
singularity near the branes, unless the combination pd(|x| 2 )/dp is well defined at p = 0. For 
the brane at p = p\ the same argument may be repeated using the coordinate p = pi—p, which 
increases as one moves away from p = 0. Keeping in mind that dp = —dp but \/—g oc +p as 
p — > 0, we see that 

- [AW A d(\ X \ 2 )/dp\^ = + [pd(| X | 2 )/dp] - =0 + [pd(|x| 2 )/dp] p=0 , (A.36) 

and so in the special case that these limits are nonsingular at the brane position (such as is 
assumed when the brane is represented as a delta-function source, or if both the background 
and perturbed geometries have only conical singularities) we see that to 2 > provided that 
(|x| 2 )' > in the near-brane limit. Furthermore, in this case u = if and only if x' = 0. 

Since the above argument assumes that the eigenmode of interest has x 7^ 0; the only 
remaining step which is required to conclude that the system is marginally stable is to sepa- 
rately show that to 2 > in the special case where x = 0- I n this case we may apply the same 
reasoning to eq. ( A.17| ), after multiplying through by _4. 3 yV 10 r*, to get 
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(A.37) 

The boundary condition can be examined as above, with ^, 3 W 10 = (.AW 4 ) 3 /W 2 oc p 3_2a °. 
Since ao = for conical backgrounds this vanishes as p 3 . We see again that u? > if 
p 3 d(|r| 2 )/dp is well-defined and non-negative at the brane position, p — > 0. In this case to = 
implies r" = T = 0. 

This leaves open only the special case of modes for which both x an d T vanish, and for 
these modes we repeat the reasoning by multiplying eq. ( |A.18| ) through by (W e /A)f*, leading 
to 
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In this case W 6 /A = (AW 4 )(W/A) 2 oc ^1+2(00-/80) x p (3-5/3 )/2_ gince p Q = 1 f or con ical 
backgrounds this varies as p" 1 as p — > 0, so if p^ 1 d(|/| 2 )/dp is well defined and non-negative 
at the brane position, p — > 0, we find uj 2 > with oj = if and only if /' = 0. 
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In this way we see how conclusions about the stability of the system can be related to 
asymptotic near-brane behaviour. Although the argument is inconclusive when the near- 
brane limit of the fluctuations is too singular, we do see that when the bulk fluctuations are 
nonsingular at the brane positions then the system is marginally stable. The stability is only 
marginal because of the known zero mode, for which T = x = while / is a nonzero constant. 



B. Comparison with previous work 



The previous stability argument, together with the one performed in comoving gauge in 
section ||, seems in contradiction with the results of ref . |14| . We therefore briefly review their 
analysis in this section and pin down the reason of the discrepancy. 

To analyze the system of linearized equations ( A.lSj , A. 17, A.16| ), ref. p4[ starts from the 
assumption that all modes are proportional to the curvature perturbation (see eq. (41) of 
Using this assumption, they find that the only existing modes are of the form 
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from which they conclude that "the spectrum of the scalar excitations consists of a zero 
mode, a first excited state with constant wavefunction and a tower of additional excited 
states with non-constant wavefunctions" . However, some other modes are present which have 



been discarded in this analysis. When plugging the relation (41) of [14] into the equations of 
motion, the authors implicitly assume that their function of proportionality cannot depend 
on the mode mass, which is an unnecessary restriction. They then derive an action (59) for 
their set of modes, but since their set of modes is not complete, the resulting action does not 
necessarily give the correct inner product for the modes. 

Even still, one could still argue that if a ghost were present among this restricted set of 
modes, it should still be present when the full set is properly considered. To see its effect, we 
concentrate for now on the ghost found in 14]. Following the prescription of [14], the ghost 
mode seems to have its origin in the two mixing modes iftx and ip2 of mass m 2 = 16g 2 which 
couple in the action as 



c *i,& = Mi {a- I6g 2 ) Vi + C h (□ - 16 5 2 ) h 
+ (□ - 16 5 2 ) + Bih (□ - IQg 2 ) 
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Due to a typo, the suggested way to diagonalize this action (ansatz (63)) does not get rid of 
the cross-terms. Instead of eq.(63), the correct sign for the change of variables should have 
been 
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with d = {C — A)/2B. This gives rise to an action where both <f) + and <p- come in with a 
positive kinetic term for any values of the parameters. To convince ourselves, we may perform 
the equivalent change of variables: 



so that the action simplifies to 



with 
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K± = AC±VACB 



— (1 + 10x 2 + x 4 ) (y(l + x 2 ) 2 + 8x 2 ± (1 + x 2 )) > 0, 

where we wrote x = ^g/q and used their expressions (60)-(62) for A,B,C, so that both JC+ 
and /C_ are positive for all values of the parameter space. Thus we do not recover the same 



result as the one argued in [ 14 ] . The presence of a ghost in their analysis seems to be due to 
a sign error in the change of variable 



In section |5.2| , we re-derive the action exactly for the entire set of modes in the theory, 
working in the more physical comoving gauge, and find the stable action which reproduces 
the field equations for the fluctuations. 



C. Decoupled Variables in Different Gauges 



In this appendix we relate the variables (x, ip, f) — in terms of which the fluctuation equations 
decouple — to the fluctuation variables in comoving and longitudinal gauges. 

Comoving gauge 

In comoving gauge we have 2A^ = and so using the gauge transformations ( |3.9[ ), we 
have 
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Longitudinal gauge 

In longitudinal gauge our notation is, $™ = —f/2, and = Vl/W — ^/2. Using the gauge 



transformation (|3.9|) then leads to 



i(0 



X 



V3_ 



2^ 



(2W 2 + X'Z')<S>® - 2U 2 ^ {1) + 2X'Z'B® 



<Po + 



2W 



+ x'z' 



4' 



2X' - 3Z' 



4V 



A y 



(0 



+ 



+ </>o 



(0 



3V 



+ 



3Z' 



2V 



(2y 
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,(2y + Z') 4V 
A a' r 



(0' 



where A^p should be related to the variables (Bm,$m,fm), (or (X(n, Tm, /m)) using the 
perturbed equation ( |A.9| ). We note that eq.( |A.9|) gives an expression for Aq up to an inte- 
gration constant f2(p). However Q, only contributes to the massless part of Ag , for which it 
can be interpreted as a redefinition of the background quantities. We can therefore set SI to 
zero in what follows. 

For ease of reference we also give here the specialization of the above relations for the 
case of conical singularities, for which we have: 

X = — X(/) 



/ 



-X 



nq \4Kq 4KqU J 



(C.2) 



_(u_ + ry\ rll) 



\2nq 2KqU J 



X' 
2nqU 



r(0' 



4W 



X 



J_ r W 

2^ 



For consistency, one can check that the field equations (3.24, 3.25, 3.26| ) in the conical case 
together with this change of variable give back the equations in longitudinal gauge (A. 16, 



A.17, A.l 



D. Some Special Functions 

In this Appendix we record some useful properties of the special functions which are used 
in the main text. In our conventions the Legendre functions are related to Hypergeometric 
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functions by 

Pu{z) = F 

Qv(z) 



-u, 1 + is; 1; -(1 - z) 



(2z)i+T (§ + 1/) 



v 1 v 3 _ 9 
1 + -,- + -;- + v.z 
2 2 2 2 



. 3 5 

u ^ , , 

T 2' 2' 



The asymptotic form of the Legendre functions as z — > ±1 can be found from well- 
established properties of the hypergeometric functions, most notably its series definition 



^ r , ab a(a + 1)6(6 + 1) z 2 



(D.l) 



provided c / 0, —1, —2, .... The potential singularities of F[a, 6; c; z\ lie at z = ±1 and z = oo, 
and our interest is in particular its behaviour at z = ±1. A useful identity for identifying 
these behaviours is 



F[a,b;c;u] = ^"^w" jjj F[a, 6; 1 + a + 6 - c; 1 - 



T(c - 6)r(c - a) 



(D.2) 



+(!-«) 



c— a— f) 



r(c)r(o + 6 - c) 



F[c — a,c — b;c — a — 6+l;l — u] 



r(a)r(6) 

This expression shows that Ffa, 6; c; z\ is finite as z — > 1 provided c > a + 6, it diverges there 
logarithmically if c = a + 6, and it can be more singular if c < a + 6: 

F[a,6;c;l] = ^ c)r(c ~ Q ~^ if c > a + 6, a,6 Z_ 
I\c — a)T(c — b) 

F[a,b;a + b;z] j-^ l og (l - z), if a, b $ Z_ 

F[-n,&;c;l] = (C ~ 6) " n £ N . 
(c)n 

Here T is Euler's Gamma function T(n) = (n — 1)!, defined for complex arguments by T(x) = 
/ °° du u x ~ 1 e~ u . (• • •)„ is the Pochhammer symbol, defined by (x) n = F p^"^ . 

For instance, using eq. (D.2) with a = —v, 6 = 1 + i/, c = 1 and u = ^(1 — z) leads 
(provided v ^ 0, 1, 2, . . .) to the following asymptotic expression 



ln[i(l + ^)] 

r(-i/)r(i + !/) 

ln[i(l + z)] 



+ O0) 



+ 0(1) as * 



(D.3) 



r(- v )T(i + v) 

If, on the other hand, v = £ = 0, 1, . . ., then the hypergeometric series terminates and Pe(z) is 
bounded at both z = ±1. This leads in the usual way to the Legendre polynomials, for which 
with Pq(z) = 1 and: 



1,2 



(D.4) 



- 34 - 



so that Pi{z) = z, P 2 (z) = ±(3z 2 - 1), P 3 (z) = \ z(5z 2 - 3), and so on. 



Similarly, using a = 1 + 6 = | + c = | + z/ and u = z 2 in the identity ( D.2 ), leads 



to 



Qu(z) 



vrr(l + u) 



( 2 ^r(i + |)r(i + |) 

± , 1 . \- In (1 - z~ 2 ) + O (1) 



In ( 1 - z" 



F 



v 1 1/ 

l + -,- + -;l;l-z~ 
2 2 2 



as z — ► ±1 , 



+ 0(1) 
(D.5) 



where we have used properties of the Gamma function to provide further simplifications. This 
function clearly has logarithmic singularities at both z = 1 and z = — 1. Notice that in the 
region of interest, —1 < z < 1, we have 1 — z~ 2 < 0, leading to a function which is complex. 
Since x ls given by the real part of Q v we have the following expressions for the near-brane 
singularities of x : 



X (z)=~hi(l-z) + 0(l) asz^l 



Ci 



C, 



r(-i/)r(i + y) 2(-i)" 



ln(l + z)+0(l) asz^-1 



(D.6) 



The homogeneous part of the equation for T is also of Hypergeometric form, leading to 
the following general homogeneous solutions: 



T h {z) = C 3 F 



V 1 V 1 n 

2 ' 2 2 ' 2' 



+ C A zF 



v 3 
2 ! 2 ; 
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Notice that F(a, 6; 6; s) = F(b, a;b;x) = 1 + ax + | a(a + l)^ 2 + • • • is independent of 6, and so 
the two Hypergeometric functions in the above expression become identical to one another in 
the special case v = 0. In fact, for v = the Hypergeometric series can be summed explicitly 
to give 



F[l,b;b;z 2 ] = Y J - 



,2/,' 



k=0 



1-Z 2 



(D. 



showing that in this case the two homogeneous solutions are (1 — z 2 ) -1 and z(l — z 2 ) -1 . 

The asymptotic behaviour of the Hypergeometric functions as z — ► ±1 is found by spe- 
cializing eq. (p.2[) to the cases (i) a = 1 + b = ^ — f,c=| and u = z 2 ; and (ii) a = | + |, 



6 = 1- 
leads to 



2' L 2 



and u = z . Keeping in mind that F[a, b ; e ; x] = 1 + abx/e + 0(x ), this 



F 
F 



V 1 ^ 1 2 

1+ 2 '2 " 2 ! 2 ;Z 

v 3 z/ 3 2 
1_ 2 '2 + 2 ! 2 ;Z 



7T 



r(i + |)r(i-|) 



i 



l-z 2 

1 



2 r(i-f)r(l + f) 



r^-iG + ln(1 ^ 2)+o(1) 

(D. 
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